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ABSTRACT 

The environment is known to affect the formation and evolution of galaxies considerably 
best visible through the well-known morphology-density relationship. It is less clear, though, 
whether the environment is equally important at a given galaxy morphology. In this paper we 
study the effect of environment on the evolution of early-type galaxies as imprinted in the 
fossil record by analysing the stellar population properties of 3,360 galaxies morphologically 
selected by visual inspection from the Sloan Digital Sky Survey in a narrow redshift range 
(0.05 ^ z ^ 0.06). The morphological selection algorithm is critical as it does not bias 
against recent star formation. We find that the distribution of ages is bimodal with a strong 
peak at old ages and a secondary peak at young ages around ~ 2.5 Gyr containing about 
10 per cent of the objects. This is analogue to 'red sequence' and 'blue cloud' identified 
in galaxy populations usually containing both early and late type galaxies. The fraction of 
the young, rejuvenated galaxies increases with both decreasing galaxy mass and decreasing 
environmental density up to about 45 per cent, which implies that the impact of environment 
increases with decreasing galaxy mass. The rejuvenated galaxies have lower a/Fc ratios than 
the average and most of them show signs of ongoing star formation through their emission 
line spectra. All objects that host AGN in their centres without star formation are part of the 
red sequence population. We confirm and statistically strengthen earlier results that luminosity 
weighted ages, metallicities, and a/Fe element ratios of the red sequence population correlate 
well with velocity dispersion and galaxy mass. Most interestingly, however, these scaling 
relations are not sensitive to environmental densities and are only driven by galaxy mass. We 
infer that early-type galaxy formation has undergone a phase transition a few billion years 
ago around z ~ 0.2. A self-regulated formation phase without environmental dependence has 
recently been superseded by a rejuvenation phase, in which the environment plays a decisive 
role possibly through galaxy mergers and interactions. 

Key words: galaxies: elliptical and lenticular, cD, galaxies: evolution, (classification, colours, 
luminosities, masses, radii, etc.), galaxies: abundances, galaxies: active, surveys 



1 INTRODUCTION 

The environment is known to be a major driver in the for- 
mation and evolution of galaxies. Its influence is best visible 
through the well-known morphology-density relationship, accord- 
ing to which early-type galaxies and morphologically undis- 
turbed galaxies are pre ferentially found in high density environ- 
ments and vice versa dOemlerl Il974l: iMelnick & Sargentl Il977l ; 
iDresslerll 1 98fA iPostman & Gellerll 1 9841) . The existence of this en- 
vironmental footprint is well estab lished, and has been repeat- 
edly confirmed both in the local iTran et all l200ll; iGoto et all 
20031 ; iHelsdon & Ponmanll2003l:lKuehn & . Rvde rj|2005t iPark et all 
20071 ; ISorrentino et al.|l2006l;lHolden et alJl200^Tvan der Wei et al.l 
20071 ; iBamford et al.l 120090 and intermediate redshift uni- 



verse 



1 Dressier et alj 1 19971; iFasano et alj bood: iTreu et al.l 2003 



I Smith et alj |2005t Ivan der Wei et al.l |2007| ; IPannella et alj I200S ; 
iTasca e t aL 2009). While the exact mechanisms responsible for this 
effect are still debated (e.g. iKodama & Smailll200ll ; iBalogh et ail 
Treu et al.l l2003), its mere existence is in good agreement 
with the idea of hie rarchical structure formation based on col d dark 
matter cosmology dBlumenthal et alj 1 9841 ; iDavis et aljl985l) . 

In contrast, it is less clear whether the environment is 
equally important at a given galaxy morphology. There is still 
major controversy about whether the formation and evolution 
of the most massive and morphologically most regular galaxies 
in the universe, i.e. earl y-type galaxie s, are affected by en- 
vironmental densities. In iThomas et al. I j2005l hereafter T05) 
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studying the local early-type galaxy population (what we call 
the 'astro-archaeology approach') we find significant evolu- 
tion as a function of environment in agreement with a large 



number of inv e stigati o ns in the literature (Trager et al. 
Longhetti et al. 2000; Poggianti et al. 2001a; Kuntschner et al. 



2000a; 



2002; Terlevich & Forbes 2002; Proctor et al. 2004; Denicolo et al 



2005 



Mendes de Oliveira et al j 



2005 



Collobertet al.l 



Schawinski et alj l2007al; ISanchez-Blazquez et al 



2006 



2006 



Rakos et alj|2007l ; lAnnibali et alj|2007l ; Ide la Rosa et alj|2007t) . In 



summary, early-type galaxies in lower densities appear younger by 
1 — 2 Gyr at given galaxy mass. 

These differences should become even more pronounced 
when p robing earlier epochs. As discussed in the review by 
iRenzinil {2006), however, studies of the intermediate redshift 
early-type galaxy population do not confirm this picture. The 
same colour evolution for bo th field and cluste r galaxies is found 
in recent redsh if t surveys jWolf et alj |2003l; iBell et al I |2004| ; 
iKoo et alj 120051 ; iTanaka et alj 120051; iBundv et al.l 120061) . This 
negligible influence from the environment is further supported 
by investigations of the fundamental plane of ellipticals, star 
formation rates, luminosity function etc. back to redshift z ~ 1, 
that again find similar evolution in clusters and in the field 



l Treu et alj 1200 ll; Ivan Dokkum & Ellis! 120031; Ivan der Wei et al] 



20051; Idi Serego Alighieri et al] 1200a |j0rgensen et al] 120061; 
van Dokkum & van der Marell 120061; di Serego Alighieri et al. 
20061; iFeulner et al] 120061; IStrazzullo et alj l2006t IPannella et ~ 



20091) . Finally, IWake et all (|2005) study clusters at 



0.3 with 

a large range of X-ray properties and also find a relatively small 
influence of the environment on the colour-magnitude relation. 

With this paper we utilise the SDSS database to check the 
result found in T05 via an increase of the sample size by more 
than an order of magnitude. The Sloan Digital Sky Survey (SDSS, 
lYork et al] 2000) provides the opportunity to explore huge homoge- 
neous samples of early-type galaxies in the nearby universe, com- 
prising several ten thousands of objects, so that a statistically mean- 
ingful investigation of the stellar population parameters of galaxies 
and their dependence on environment can be attempted. We analyse 
the stellar population parameters luminosity-weighted age, metal- 
licity, and a/Fe element ratio of 3,360 early-type galaxies drawn 
from the SDSS in a narrow redshift range (0.05 ^ z sC 0.06). 
A major strength of our approach is the application of strict mor- 
phological selection criteria based on visual inspection of 48,023 
galaxies. We confirm the basic result of T05 as far as 'archaeo- 
logical downsizing' and the major formation epochs of early-type 
galaxies are concerned. However, different from T05, the latter turn 
out to be insensitive to changes of the environment for the bulk of 
the early-type galaxy population. 

The paper is organised as follows. In Section 2 we describe in 
brief the construction of the galaxy catalogue, and the method to 
derive stellar population parameters. The results on stellar popula- 
tion scaling relations are presented in Section 3, galaxy formation 
epochs are discussed in Section 4. The paper concludes with dis- 
cussion of the results and literature comparison in Sections 5 and 6. 



2 DATA SAMPLE AND ANALYSIS 

The sample utilised here is part of a project called MOSES: 
Morphologically Selected Early-types in SDSS. We have collected 
a magnitude limited sample of 48,023 galaxies in the redshift range 
0.05 ^ z ^ 0.1 with apparent r- band magnitude brighter than 16. 8 
from the SDSS Data Release 4 jAdelman-McCarthv et aT]|2006l) . 



Model magnitudes are used throughout this work. The most radical 
difference with respect to other galaxy samples constructed from 
SDSS is our choice of purely morphological selection of galaxy 
type. All galaxies in this sample have been visually inspected and 
classified by eye into early- a nd late-type morphology. We have al- 
ready used this approach in ISchawinski et al] ]2007bl) . and more 
details on the morphological selection and other sample charac- 
teristics can be found there. In the following, we provide a brief 
summary with focus on the aspects that are most relevant for the 
present work. 



2.1 Morphological selection 

We visually inspected the 48,023 objects and divided the sam- 
ple in 31,521 late-type (spiral arms, clear disc-like structures) and 
16,502 early-type (roundish, elliptically shaped) galaxies. Hence, 
34 per cent of the objects are classified as 'early-type'. This is in 
good agreement with the typical value found in average environ- 
ments (see references in the Introduction). For the present study 
we aim to minimise contamination by late-type galaxies, hence our 
relatively strict selection criteria yield a conservative early-type 
galaxy fraction that is slightly smaller th an the recent determina- 
tion based on SDSS bv lvan der Wei et alj (~ 40 per cent. 120071) but 
in good agreement w ith other determinations in the literature (see 
iBamford et alj2009l and references therein). 

Elliptically shaped objects with irregular structures, poten- 
tially indicative of tidal tails originating from recent galaxy-galaxy 
interaction/merger events, were included in the 'early-type' cate- 
gory. The principal aim was to separate spiral and disc-like galaxies 
from early-types. The major advantage of this strategy is that our 
sample is not biased against star forming elliptical galaxies. No fur- 
ther selection criteria based on stellar population properties such as 
colour, spectral shape, or concentration index have been applied. 
Indeed, this catalogue differs mos t from previous S PSS-based 
early-type galaxy catalogues (e.g.. iBernardi et al]|2003h through 
its inclusion of (typically low-mass) objects with unusually blue 
colours. We have ex tended this approa ch recently through the 
Galaxy Zoo project dLintott et alj 120081) . which uses several ten 
thousands of volunteers on the internet to classify galaxies through 
a web interface. Among other projects, these classifications have 
led to the construction of a large sample of blue early-type galaxies 
(Schawinski et al 2009b; see also Kannappan et al 2009). 



2.2 Environmental density 

For the estimate of the environmental densit y we adopt the method 
dev eloped and describ ed in full detail in lSchawinski et al] d2007al) 
and Yoon et al. 1 120081) . In a nutshell, we calculate the local number 
density of objects brighter than a certain absolute magnitude in a 
6 Mpc sphere around the object of interest. A Gaussian weight 
function centred on the object ensures to avoid contamination from 
neighbouring structures, and an adjustment of the side length of the 
'sphere' in redshift direction is invoked in order to account for the 
'finger of God effect' . We have tested various sizes for the Gaussian 
weight functions and found that the final results of this work are not 
sensitive to the particular size of the sphere around 6 Mpc. Note 
that the results using spheres of only a few Mpc size were noisier 
because of the sparse sampling by SDSS fibers. The 6 Mpc chosen 
here can be considered conservative for the purposes of this work in 
that this choice assures the definition of low environmental density 
to indicate truly isolated regions. 
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We obtain the 3-dimensional volume density at the location 
of each object, hence an estimate of the local environmental den- 
sity. This process does not constrain the overall structure in which 
the object resides, but is a sophisticated number density providing 
a measure of the number and proximity of galaxies around a point 
in space. In other words, we cannot distinguish the outskirts of a 
galaxy cluster from a group environment of the same local density. 
We decided not to include further c onstraints like cross-c orrelation 
with cluster catalogues like in, e.g jBer nardi et al] J200(J) . because 
also this method suffers from incompleteness problems. More gen- 
erally speaking, we did not attempt to correlate our density esti- 
mates with physical structures, as we are not aiming for classifica- 
tions such as ' field' , group' , or ' cluster' . In any case, as discussed in 
Section 5, our results are overall consistent with other SDSS studies 
in the literature as far as the general influence of the environment is 
concerned. 



2.4 Stellar population parameters 



2.3 Kinematics and Lick line index measurements 

The kinematics of gas and stars are determined for al l 48,023 ob - 
jects, using the code GANDALF developed bv lSarzi et al.l J2006I) . 
Stellar population and emission line templates are fitted simultane- 
ously to the galaxy spectrum. The outcome are emission line fluxes 
and kinematics of the gas. On the stellar population side, we obtain 
li ne-of-sight velocity dispersions based on the pixel fitting method 
of ICappellari & Emselleml j2004) with a typical error of 0.023 dex 
in log o . The error distribution deviates from a Gaussian with a 
sharp cut-off at 0.01 dex and a tail extending to larger errors of 
~ 0.1 dex. Through the subtraction of the emission line spectrum 
from the observed one, we get the clean absorption line spectrum 
free from emission line contamination. This spectrum is used for 
any further analysis of absorption features. The stellar template 
spectrum is not used for scientific analyses. 

On each spectru m, we measure the 25 standard Lick ab- 
sorpti on line indices dWorfhev et all 1 19941 ; | Worthev & Ottaviani 
1997) following the most recent index definitions of iTrager et al. 
(1998). For this purpose the spectral reso lution is reduced to 
the w avelength-dependent Lick resolution dWorfhev & Ottavianil 
1 19971) . The measurements are then corrected for velocity disper- 
sion broadening. The correction factor is evaluated using the best 
fitting stellar template and velocity dispersion obtained previously. 
Errors are determined by Monte Carlo simulations on each spec- 
trum individually based on the signal-to-noise ratios provided by 
the SDSS. 

It should be noted that an important step in the measurement 
of Lick indices cannot be undertaken. Ideally, Lick standard stars 
need to be observed with the same instrument set-up, in order to de- 
termine (and correct for) offsets caused by the shape of the spectral 
energy distribution. This is in principle necessary, as the Lick/IDS 
system is not based on flux-calibrated spectra. Such calibration 
stars at sufficiently high signal-to-noise are not available in SDSS. 
However, it ought to be expected that the shape of the spectrum is 
of little importance for absorption indices where line feature and 
pseudo-continuum windows are narrow and close enough in wave- 
length. This condition is actually fulfilled for most Lick indices 
with the pseudo-continuum bands being placed directly adjacent 
to the line feature itself. A detailed analysis of Lick standard stars 
extracted from the SDSS archive sugg ests that devia tions from the 
true Lick system are indeed negligible iCarsonll2007h . 



We use the stellar population models of lThomas et alj d2003l . [2004L 
hereafter TMB) to derive luminosity-weighted ages, metallicities, 
and a/Fe ratios from the Lick absorption lin e indices by means 
of the minimised x 2 technique introduced by IProctor & Sansoml 
d2002l) . The main motivation for this choice is that the \ 2 fitting is 
very efficient and fast, and can easily be applied to large data sets 
like SDSS. A further advantage is that it considers a large wave- 
length range of the spectrum and hence gives a better estimate of 
the real average population. Finally, the consideration of a larger set 
of absorption features increases the signal-to-noise ratio. Errors are 
obtained from the probability function based on the x 2 distribution. 
In practice, we search for the combination of those three parameters 
that minimise the \ 2 between model and observational data for all 
25 indices observed. This is different from the approach followed 
in T05, where only a selection of well calibrated indices is consid- 
ered to constrain the stellar population parameters directly rather 
than in a statistical fashion. It has been shown, however, that these 
two o rthogonal methods yield consistent results dThomas & Daviesl 
l2006h . 

We follow an iterative approach. In the first step we fit all 25 
indices. We compute the probability distribution for different val- 
ues of x 2 at its minimum using the incomplete T function. This 
yields the probability P that the observed x 2 f° r a correct model 
should be less than the value x' 2 obtained in the fit. If P ^ 0.999 
we deem the fit unacceptable. In this case we discard the index 
with the largest x 2 from the fitting procedure, re-do the fit, and 
re-calculate P. This procedure is iterated until P < 0.999. The 
indices Ca4455 and NaD were discarded in more than 10 per cent 
of the cases, and were hence generally excluded. Their exclusion 
is furthe r supported by the fact that these indices a re not well cal- 
ibrated (Maraston 1-t~al1l2003l : iThomas et alj|2003l) . In more than 
half of the cases (65 per cent), no further index was removed (hence 
all 23 indices were used). In 20, 8, and 4 per cent of the objects 1, 
2, and 3 indices were removed, respectively. In other words, for 97 
per cent of all objects, at least 20 indices were included in the fit. 

The average errors in log(age), metallicity, and [a/Fe] ra- 
tio are 0.10 dex, 0.07 dex, and 0.06 dex, respectively. The er- 
ror distributions are well approximated by Gaussian functions with 
standard deviations of 0.04, 0.02, and 0.03 dex, respectively. Fi- 
nally it should be noted that size of the SDSS fibers used for the 
spectroscopy is relatively large (3 arcsec). The stellar population 
properties discussed here therefore refer to the global stellar pop- 
ulations in the galaxies roughly within or only slightly below the 
half-light radius. As early-type galaxie s are known to have metal- 
licity gradients (e.g. iDavies et alj|l993h the effect of these when 
using light-average quantities is that we slightly underestimate the 
true age of the stellar population, simply because the metal poor 
halo provides blue light which is interpreted as a slightly younger 
population. We expect this effect to be small, though, because the 
amount of metal-poor stars in luminous re d galaxies in SDSS i s 
of the order of only a few per cent by mass lMara stonetaiT2 009). 
Moreover, the a/Fe ratio is not affected by this as early-type galax- 
ies e xhibit no significant gradient in a/Fe within the effective ra - 
dius dMehlert et aLll2003l ; lAnnibali et alj|2007l iRawle et alj|2008h . 
Note also that no gradient in ag e is typically found in local early- 
type galaxies dSaglia et alJl200d:lTamura et alJl2000l:lMehlert et all 
2003 ; IWu et alfeOOShlAnnibali et alj2007l ; lRawle et alj|2008h . 
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Figure 1. Left-hand panel: Distribution of velocity dispersion for different sample selections in redshift. The grey shaded histogram is the entire sample of 
16,502 early-type galaxies in the redshift range 0.05 ^ z ^ 0.1. Blue and red histograms show low and high redshift selections, respectively. Right-hand 
panel: Contours of light-averaged age as function of velocity dispersion. Colours as in left-hand panel. The final sample in this study is selected from the 
lowest redshift range (0.05 ^ z ^ 0.06, blue lines). 



2.5 Final selection 



3.1 Previous work 



We applied further selection criteria to the 16,502 early-type galax- 
ies in our catalogue. The redshift range sampled (0.05 ^ z ^ 0.1) 
is small and corresponds to a time window of only ~ 600 Myr. 
Still, this age difference produces significant selection effects. The 
left-hand panel of Fig. [T] shows the distribution of velocity dis- 
persions for low redshift (0.05 ^ z ^ 0.06, blue line) and high 
redshift (0.09 ^ z ^ 0.1, red line) selections in comparison 
with the entire sample (grey shaded histogram). It can be seen 
that despite the small range in redshift, completeness problems at 
the low-u end get worse with the inclusion of higher redshift ob- 
jects. This leads to a very bad sampling of the age-cr relationship, 
as shown by the right-hand panel of Fig. [T] On top of this selec- 
tion effect in a, small evolutionary effects complicate the deriva- 
tion of stellar population parameter s in spite of the small time win- 
dow sampled dBemardi et alj|200of) . We therefore took a conserva- 
tive approach and decided to consider only the lowest redshift bin 
(0.05 ^ z ^ 0.06, blue lines in Fig. [TJ providing an acceptable 
coverage in velocity dispersion down to logcr/km/s ~ 1.9. This 
results in a total of 3,360 early-type galaxies. This redshift range 
corresponds to a time interval of only 100 Myr at an average look- 
back time of 700 Myr adopting f2 m = 0.24, fl A = 0-76, and 
H = 73 km/s/Mpc jTegmark et al.ll2006] : |per"cival et all2007h . 



3 STELLAR POPULATION SCALING RELATIONS 

The characteristic mass of a galaxy populat ion increases with en- 
vironmental density (e.g., .Baldry et al., 2006), indicating that lower 
density environments host less massive galaxies. This mass bias 
has to be eliminated for a meaningful study of the influence of en- 
vironment. We therefore investigate the environmental dependence 
of the correlations of the stellar population properties with velocity 
dispersion a and dynamical galaxy mass Md yn . From these 'stellar 
population scaling relations' we derive formation histories follow- 
ing the method of T05. 



While work on the classical scaling relations of early-type galaxies 
such as colour-magnitude relation, Fundamental Plane, and Mg- 
a relation has suggested relatively old ages and high for mation 
redshifts already almost two decades ago (see review by iRenzinil 
2006, and references therein), the efforts of the last decade to find 
correlations between galaxy mass and the three stellar population 
parameters age, metallicity, and a/Fe ratio have succeeded only 
recently. M ostly because o f the well-known age-metallicity degen- 
eracy (e.g. IWorthevll 19941) , the direct determination of luminosity- 
weighted ages turned out to be extremely difficult. When ages 
of early-type galaxies were explicitly derived, no evidence for 
a corre lation with mass , luminosity or velo c ity dispersion was 



a correlation with mass, luminosity or velocity dispersion was 
found Jj0rgensenl [l99p ; iTrager et alj 12000a 1 Kuntschner l20ori 



Kuntschner et al.1 l200lt IPoggianti et alj 1200 lbM IVazdekis et all 



200 it iTerlevich & Forbesll2002h IMehlert et aUl2003l). As the most 
notice able exception, the study of ICaldwell. Rose & Concannonl 
(2003) hints for the presence of an age-sigma relationship, but 
driven mainly by the lowest-mass objects and with extremely large 
scatter. 

Somewhat more convincing correlations with velocity disper- 
sion were found for metallicity and the a/Fe element abundance 
ratio in the studies mentio ned above. In particular the latter could 
be further improved (T05 ; ITrager et alj|20003 ; IProctor & Sanson] 
l2002t IMehlert et al] [2003) after element abundance ratio sensitive 
stellar population models were used (TMB: ITrager et al.l l2000b). 
Besides allowing for a quantitative determination of element ratios, 
the inclusion of element ratio effects allows for a clean definition 
of total metallicity and in this way helps improving the accuracy of 
the age determination. 

A major obstacle, however, is that most indices depend on all 
three stellar population parameters. This leads to tight correlations 
between the er rors of these parameters, which hampers the ir deter- 
mination (T05: lTrager et al.l2000iJ ; lKuntschner et al.l200lh . In par- 
ticular small data sets are severely affected by this problem. T05 
follow therefore a Monte Carlo approach and seek correlations be- 
tween stellar population parameters and galaxy mass through com- 
parison of mock galaxy samples with observational data. Through 
this method, well-defined relationships between all three stellar 
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Figure 2. Contour plots of the relationship between stellar velocity dispersion and luminosity-weighted age for various environmental densities as indicated 
by the inset histograms. The environmental density is proportional to the number density per volume, but no precise physical units are associated to it. The 
lowest and highest density bins contain 658 and 571 out of 3,360 galaxies, respectively. The dashed line separates an old red sequence population (orange 
contours) from rejuvenated objects in the blue cloud with light-averaged ages smaller than 2.5 Gyr (cyan contours). The fraction of this latter population is 
given by the label. The underlying grey contours include both populations. The solid line is a linear fit to the red sequence population, the parameters of the fit 
are given at the top of each panel. The dotted line is the fit for all environmental densities (parameters from top left-hand panel). Its distribution is shown by 
the top right-hand panel (same colour coding). The label gives the standard deviation a for the fit. The age-cr relationship for the red sequence population is 
independent of environment, while the rejuvenation fraction increases with decreasing density. 



population parameters (including luminosity-weighted age) and 
galaxy velocity dispersion could be identified. Other recent work 
on extremely large data samples (mostly SPSS) basically confirm 
these results (Nelan et al. 2005; Bernardi et al. 2006; Gallazzi et al. 
2006; Clemens et al. 2006; Jimenez et al. 2007; G raves et al . 2007; 
Allanson et all2009h . 

3.2 Age as a function of environment 

Fig. [2] shows contour plots for the relationship between stellar ve- 
locity dispersion and luminosity-weighted age for various envi- 
ronmental densities as indicated by the inset histograms. The top 
left-hand panel includes all environments. The three bottom panels 
show three different choices of environments (orange histograms), 
i.e. the median density around the peak of the distribution (middle), 
and the lowest (left-hand panel) and highest (right-hand panel) den- 
sity ends of the distribution. The latter are expected to show truly 
isolated and densest cluster galaxies and contain 658 and 571 out 
of 3,360 galaxies, respectively. Note that our density estimator is 
proportional to the number density per volume but does not have 
a strictly physical meaning, hence we do not attempt a more accu- 
rate identification with astrophysical objects like galaxy groups and 
clusters. 

The top-left and right-hand panels show that the distribution of 
ages is bimodal with a major peak at old ages (orange contours) and 



a secondary peak at young ages around ~ 2 Gyr (cyan contours) in 
analogy to 'red sequence' and 'blue' cloud identified in galaxy pop- 
ulations (see also Fig. [5}. The sub-population of young early-type 
galaxies could be identified here mostly thanks to the purely mor - 
phological selection algorithm (see also lSchawinski et al. 2009b). 
The dashed line at log (age) = 0.4 divides the two populations. 
The objects below this line have very young light-averaged ages 
(< 2.5 Gyr), and have most likely been rejuvenated by recent mi- 
nor star formation events at lookback times < 2.5 Gyr. Note that 
a minor rejuvenation event adding only a few per cent of young 
stella r populations has a large effect on the luminosity-average age 
(e.g.. lGreggiolll997l ; [Serra & Tragej|2007h . Therefore, we call the 
objects below the dashed line 'rejuvenated', and the fraction of ob- 
jects in this category 'rejuvenation fraction'. The latter is about 
10 per cent as indicated by the label in the figure. It can also be 
seen that rejuvenation occurs mostly in lower-mass galaxi es with 
log o < 2.2 as previously found in other studi es (T05, lYi et al.l 
l2005l ; lThomas & Daviesll2006l ; lDonas et all2007l) . 

We discuss the rejuvenated population in the 'blue cloud' sep- 
arately and derive stellar population scaling relations for the 'red 
sequence' population that dominates the mass budget of the galaxy 
population. The linear fit to this population is shown by the solid 
line. The parameters of the fit are given at the top of the panel. 
The distribution about this fit is shown in the top right-hand panel 
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Figure 3. Contour plots of the relationship between stellar velocity dispersion and a/Fe ratio for various environmental densities as indicated by the inset 
histograms. The environmental density is proportional to the number density per volume, but no precise physical units are associated to it. The lowest and 
highest density bins contain 658 and 571 out of 3,360 galaxies, respectively. The cyan contours are rejuvenated objects with light-averaged ages smaller than 
2.5 Gyr (see Fig. [2), the orange contours are the old red sequence population. Underlying grey contours include both populations. The solid line is a linear 
fit to the red population, the parameters of the fit are given at the top of each panel. The dotted line is the fit all environmental densities (parameters from top 
left-hand panel). Its distribution is shown by the top right-hand panel (same colour coding). The label gives the standard deviation a for the fit. The a/Fe-a 
relationship for the red sequence population is independent of environment. 



(black histogram and orange line). The cyan line is the rejuvenated 
population and is clearly offset. 

The fits to the red sequence population are repeated in every 
environment bin. The dotted lines in all panels indicate the fit in- 
cluding all environmental densities. It can be seen clearly from the 
plots that the resulting fit parameters are consistent with no varia- 
tion as a function of environment within their 1-a error bars. We 
conclude there is no considerable change as a function of environ- 
mental density, hence the age-a relationship is independent of en- 
vironment. 

Very different is the behaviour of the blue cloud (rejuvenated) 
population. In this case, the environment plays a critical role. The 
rejuvenation fraction increases with decreasing environmental den- 
sity. Hence early-type galaxies in lower density environments are 
not generally younger (at a given mass), but the fraction of re- 
juvenated galaxies is higher. This depe ndence on environment is 
very similar to what has been found bv lSchawinski et al] |2007a) 
albeit tracing only the most recent residual star formation episodes 
in early-type galaxies through near-UV colours. 

Blue horizontal br anches from metal-poor or peculiar metal- 
rich s ub-popul a tions dGreggio & Renzi nf Il99d : | Dorman et all 
1 19951 : I Yi et al J Il998|) may mimic an apparent presence of 
young populations {M araston & Thomas! l2000l : iLee et all |2000| : 
iMaraston et"ai] 12001 iTrager et all 120051 T05). This degeneracy 
is very difficult to disentangle in unresolved stellar populations. 
We use therefore further indicators of residual star formation that 



are independent of horizontal branch morphology for consistency 
checks. The a/Fe element abundance ratio is a useful quantity, be- 
cause it quantifies the relative importance of delayed chemical en- 
richment fro m Type la supernovae owing to late residual star for- 
mation (e.g., Greggio & Renzini 1983; Matteucci & Greg giol 19861 : 
lMatteuccilTi994 lThomasetalJll999l : (Greggioll2005h . We further 
study emission line properties as indicators for possible star for- 
mation and/or AGN activity, as the ancillary use of emission line 
diagnostics allows us to break this degeneracy. 

3.3 The a /Fe ratio as a function of environment 

Fig. [3] shows contour plots for the a/Fe-a relations using the same 
symbols and colour coding as in Fig. [2] We confirm previous find- 
ings of a very well-defined and tight correlation between veloc- 
ity dispersion and a/Fe ratio. The linear fits are again performed 
on the red sequence population (orange contours) as defined in the 
previous section. The distribution about this fit is shown in the top 
right-hand panel (black histogram and orange line). The cyan line 
is the rejuvenated population. The dotted lines in all panels indi- 
cate the fit including all environmental densities. Also in this case 
the slope and zero-point of the correlation are consistent with no 
variation as a function of environment. 

The environment-dependent rejuvenated population (cyan 
contours) is offset from the general fit to slightly lower [a/Fe] ra- 
tios by ~ 0.1 dex. Note that errors in age and a/Fe do not correlate 
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Figure 4. Contour plots of the relationship between stellar velocity dispersion and total metallicity for various environmental densities as indicated by the inset 
histograms. The environmental density is proportional to the number density per volume, but no precise physical units are associated to it. The lowest and 
highest density bins contain 658 and 571 out of 3,360 galaxies, respectively. The cyan contours are rejuvenated objects with light-averaged ages smaller than 
2.5 Gyr (see Fig. [2), the orange contours are the old red sequence population. Underlying grey contours include both populations. The solid line is a linear fit 
to the red population, the parameters of the fit are given at the top of each panel. The dotted line is the fit for all environmental densities (parameters from top 
left-hand panel). Its distribution is shown by the top right-hand panel (same colour coding). The label gives the standard deviation <r for the fit. The Z/H-cr 
relationship for the red sequence population is independent of environment. 



such that they could create the offset in a/Fe seen here (see Fig. 
3 in Thomas et al 2005). This result is critical, as it supports the 
interpretation that the low light-averaged ages of the rejuvenated 
population are indeed caused by the presence of residual star for- 
mation in these objects. This conclusion is further supported by the 
fact that rejuvenated elliptical galaxies with low a/Fe ratios tend to 
have extra lig ht in the centre pointing toward a dissipational forma- 
tion scenario dKormendv et alj|2009t) . Fig. [3] shows that the offset 
in a/Fe ratio is relatively small, however, which suggests that only 
minor star formation episodes are responsible for this rejuvenation. 
This is in agreement with the mean stellar ages derived from the 
mid-infrared fluxes of early-type galaxies that onl y allow for small 
amounts of recent star formation dTemi et alj2005h . 



Like for age and a/Fe ratio, these relationships are independent of 
environment. 

Fig. E] further shows that the rejuvenated galaxies (cyan sym- 
bols) have slightly higher metallicities than the bulk population by 
~ 0.1 dex. This agrees well with previous findings of an age- 
metallicity anti-corr elation at given a in early-type galaxies (e.g., 
iTrager et al.l 12000 j) . The slight increase in metallically suggests 
that the residual star formation in these galaxies does not involve 
purely metal-poor pristine gas but must include at least some frac- 
tion of previously enriched interstellar medium, either from internal 
re-processing or external accretion. Again, the fact that the offset is 
small suggests that the rejuvenation event involves only a minor 
fraction in mass. 



3.4 Total metallicity as a function of environment 

The contour plots for the Z/H-a relations are shown in Fig. fus- 
ing the same symbols and colour coding as in Fig. [2] Again, we 
confirm the existence of a well-defined Z/FL-o relation. Our data 
mainly improve on its tightness and significance. The linear fits are 
again performed on the red sequence population (orange contours) 
as defined in the previous section. The distribution about this fit 
is shown in the top right-hand panel (black histogram and orange 
line). The cyan line is the rejuvenated population. The dotted lines 
in all panels indicate the fit including all environmental densities. 



3.5 Colour 

As a cross-check of the stellar population parameters derived from 
absorption line indices we briefly discuss u—r colour. Fig.|5]shows 
contour plots for the u — r-a relation using the same symbols and 
colour coding as in Fig. [2] All environmental densities are included. 
The linear fits are again performed on the red sequence population 
(orange contours) as defined in the previous section. The rejuve- 
nated population (cyan contours) with light- averaged ages below 
2.5 Gyr as defined in Fg.|2]is well separated from the red popu- 
lation showing bluer it — r colour by about 0.2 mag. The old pop- 
ulation forms a well-defined red sequence, while the rejuvenated 
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Figure 5. Contour plots of the relationship between stellar velocity disper- 
sion and u—r colour including all environmental densities. The cyan con- 
tours are rejuvenated objects with light-averaged ages smaller than 2.5 Gyr 
(see Fig. 0, the orange contours are the old red sequence population. Un- 
derlying grey contours include both populations. The solid line is a linear 
fit to the red population, the parameters of the fit are given at the top of the 
panel. The rejuvenated population is clearly visible through its blue u — r 
colour, which is offset from the general trend by ~ 0.2 mag. 



objects occupy the region that is known as the 'blue cloud'. The 
u—r colour clearly reinforces the conclusions drawn from Figs. [2] 

3.6 Star formation activity and AGN 

It is interesting to analyse the ionisation states of the gas in the 
early-type galaxies as a function of their position on the age-cr re- 
lationship. In particular, the possibility of blue horizontal branch 
stars mimicking the presence of young stellar populations can be 
ruled out if rejuvenation coincides with star formation activity 
detected through emission lines. The majority of the sample (77 
per cent) shows no significant emission lines, hence are classified 
as 'passive'. To the remaining 'active' 23 per cent we apply the 
BPT classification scheme that separates star format i on from AGN 
jBaldwin et alj Il98ll: IVeilleux & Osterbrockl 1 19871 : iKewlev et all 
l200ll ; iKauffmann et alj|2003k iMiller et alj|2003l). A detailed BPT 
analysi s of the present sample is presented in Schawinski et al. 
(2007b) and we refer the read er to this paper for furth er details. 
Using the classifications from Schawinski et "all i l2007bl) . we find 
that 16 per cent of the objects with emission lines are star forming, 
21 per cent are star forming/AGN composites, 16 per cent have 
AGN Seyfert-like emission lines, and 47 per cent show LINER- 
type emission. 

In Fig. [6] we plot contours for these sub-populations over the 
age-cr relationship of the whole sample (analogue to Fig. [2}- Star 
forming objects are plotted in blue, star forming/AGN compos- 
ite objects in green, AGN with Seyfert-like emission in orange, 
and finally AGN with LINER-like emission and passive galaxies 
are plotted in red contours. It turns out that all objects that show 
signs of star formation activity (blue and green contours) have re- 
juvenated stellar populations. Some fraction of these are star form- 
ing/AGN composites, but all objects that host AGN in their centres 
without any star formation harbour old stellar populations. More 
specifically, we have detected star formation activity in 286 out of 
3360 objects (8.51 per cent) compared to 341 (10.15 per cent) that 
have been classified as rejuvenated. This implies that the majority 



Figure 6. Contour plots of the relationship between stellar velocity dis- 
persion and age including all environmental densities (analogue to Fig. [2). 
The coloured contours are sub-populations separated through their emis- 
sion line classifications star forming (blue), star forming/AGN composite 
(green), AGN Seyfert-like emission (orange), AGN LINER-like emission 
and passive (red). The underlying grey contours include all galaxies. All 
objects with traces of ongoing star formation (blue and green contours) be- 
long to the rejuvenated population defined in Fig. [2] while galaxies with 
AGN have considerably older stellar populations. 



of the rejuvenated galaxies (84 per cent) are actively forming stars, 
which is well consistent with the detection of significant amount s 
of neutral hydrogen in low-mass ellipticals iMor ganti et alj|2006f) . 
This clearly rules out the possibility of blue horizontal branch stars 
mimicking the presence of young stellar populations. 



3.7 Dynamical masses 

So far in this paper we have used line-of-sight stellar velocity dis- 
persion a as a proxy for galaxy mass. The major reason is that 
cr is independent of the stellar population properties derived here 
and can be easily measured in gal axy spectra from the broaden- 
ing of absorption line features (e.g jFaber & Iacksonll976llBender] 
1990). However, a is a difficult quantity to predict in galaxy for- 
mation models, the direct output of which are stellar and dynami- 
cal masses. Hence, we determine dynamical masses from the scal- 
ing between virial dynamical galaxy mass Md yn with stellar veloc- 
ity dis persion a, and effective radius provided by Cappellari et al. 
(2006). Effective radii are taken from the SDSS database. The 
typical error in log Md yn is 0.05 dex with a tail extending to 
~ 0.2 dex. Note that this 'dynamical' mass should be a good ap- 
proximation also of the baryonic mass, as baryonic matter is dom- 
inating ttejnass_budgetjnearlyj^y^ the effective 
radius dCappellari et alj200d : lThomas et alj2007l) . This is different 
from the approach in T05, where a scaling with observed luminos- 
ity through model M/L ratios as a function of age and metallicity 
is used. We prefer the dynamical mass estimates, as they are inde- 
pendent of the stellar population parameters derived. 

Fig.|7]presents the equivalent plots to Figs.[2][3] and|4]showing 
the resulting relationships between dynamical mass and stellar pop- 
ulation parameters light-averaged age (left-hand panel), a/Fe ele- 
ment ratio (middle panel), and total metallicity (right-hand panel). 
All environmental densities are included. The dashed line (left- 
hand panel only) separates again the red sequence population (or- 
ange contours) from the rejuvenated objects with light-averaged 
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Figure 7. Contour plots of the relationship between dynamical mass and light-averaged age (left-hand panel), a/Fe element ratio (middle panel), and total 
metallicity (right-hand panel) including all environmental densities. The dashed line (left-hand panel) separates the old red sequence population (orange 
contours) from the rejuvenated objects with light-averaged ages smaller than 2.5 Gyr (cyan contours). The fraction of this latter population is given by the 
label. The underlying grey contours include both populations. The solid line is a linear fit to the red population, the param eters of the fit are given at the top of 
each panel. The dynamical mass is calculated from the scaling with velocity dispersion and effective radius adopted from Cappellari e t al.l j200d) . 



ages smaller than 2.5 Gyr (cyan contours), while the underlying 
grey contours include both populations. The solid line is a linear fit 
to the red sequence population, the parameters of this fit are given 
at the top of each panel. It is evident that the scaling relations are 
just as well defined for dynamical mass as they are for velocity 
dispersion. 



3.8 Rejuvenation fractions 

The results presented in this work show that it is the fraction of re- 
juvenated galaxies (galaxies in the blue cloud) that is dependent 
of environment, while the stellar population scaling relations of 
the red sequence galaxies are environment independent. So far, we 
have not disentangled this from the general correlation between age 
and mass. As the rejuvenation fraction increases also with decreas- 
ing velocity dispersion and galaxy mass, its environmental depen- 
dence may just be caused by mass segregation, which is the fact that 
more massive galaxies live preferentially in denser environments. 

Fig. [8] quantifies this in more detail. The plots show the frac- 
tion of rejuvenated galaxies as a function of environmental density 
(left-hand panel) and velocity dispersion (or galaxy mass, right- 
hand panel) for three different bins in velocity dispersion and en- 
vironmental density, respectively. Fig. [^demonstrates that the frac- 
tion of rejuvenated early-type galaxie s increases as a functio n of 
both decreasing galaxy mass (see also Schawinski et al. 2006) and 
decreasing environmental density dSchawinski et alj |2007a). with 
the effect of environment being slightly weaker. Fractions vary 
from only a few per cent for the most massive galaxies in dens- 
est regions up to 45 per cent for the smallest ellipticals (a few 
10 10 Mq) in the least dense regions of the universe. It can be 
seen that the rejuvenation fractions as well as their dependence on 
environment are largest for the smallest galaxies (blue line in the 
left-hand panel). The most massive early-type galaxies have experi- 
enced no rejuvenation event independent of environment (magenta 
line in the left-hand panel). The rejuvenation fraction among low- 
and intermediate mass galaxies (blue and orange lines) increases 
from the highest to the lowest density bin at > 2a significance. 
The increase of the rejuvenation fraction as a function of velocity 
dispersion is most pronounced for the lowest environmental densi- 
ties (blue line in right-hand panel). 



3.9 The final scaling relations 

The age fits for the red sequence population shown in Figs. [2] and 
[7] as orange solid lines are log age(Gyr) = —(0.25 ± 0.05) + 
(0.52 ± 0.02) log a and log age(Gyr) = -(0.76 ± 0.09) + 
(0.15 ± 0.01) log Mdyn- Note, however, that the sample discussed 
here covers the redshift range 0.05 ^ z ^ 0.06. Even though 
these are very moderate redshifts, they correspond to look-back 
times that are sig nificant for the purposes of this work (see also 
lBernardietai]|2006l) Adopting Qrn = 0-24, fl A = 0.76 and 
H = 73 km/s/Mpc dTegmark et alJl200d : |Percival et alj|2007l) we 
obtain an average lookback time of 700 Myr for this redshift inter- 
val. By correcting for this time evolution (on the linear age scale), 
the zero points of the above relationships (on logarithmic scale) 
increase to —0.11 and —0.53, and the slopes decrease to 0.47 and 
0.13. We assume that both metallicity and a/Fe ratio do not evolve. 

We then obtain the following local stellar population scaling 
relations, age (log (t/Gyr)), total metallicity (Z/Fi), and a/Fe ele- 
ment ratio as functions of velocity dispersion a (in km/s) and dy- 
namical mass Afdyn (in Mq). Note that these are independent of 
the environment. 



log t = 


-(0.11 ± 


0.05) 


+ 


(0.47 ± 0.02) 


log(o-) 


[Z/H] = 


-(1.34 ± 


0.04) 


+ 


(0.65 ±0.02) 


log(cr) 


a/Fe] = 


-(0.55 ± 


0.02) 


+ 


(0.33 ±0.01) 


log(a) 


log t = 


-(0.53 ± 


0.09) 


+ 


(0.13 ±0.01) 


log(M) 


[Z/H] = 


-(2.40 ± 


0.07) 


+ 


(0.22 ±0.01) 


log(M) 


a/Fe] = 


-(0.95 ± 


0.04) 


+ 


(0.10 ±0.01) 


log(M) 



These relationships represent the red sequence population and 
are valid for a large range in galaxy masses (o > 100 kms -1 or 
Mdyn = 3 • 10 10 Mq) and all environmental densities. Note that 
a major fraction of galaxies deviates from these scaling relations at 
the low ends of the mass and environmental density distributions, 
as the fraction of rejuvenated galaxies increases to ~ 45 per cent 
(see Fig.[8]>. 

The observed spread in log age, metallicity, and [a/Fe] ratio 
in Eqn.Qjis on average 0.23 dex, 0.11 dex, and 0.07 dex, respec- 
tively (see top right-hand panels in Figs. [2] [3] and|4]l. Subtracting 
the errors estimated in Section [2] we obtain an average intrinsic 
scatter of 0.21 dex, 0.08 dex, and 0.02 dex, respectively. The in- 
trinsic scatter in age increases slightly to 0.27 dex at the lowest 
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Figure 8. Left-hand panel: Fraction of rejuvenated galaxies as a function of environmental density (left-hand panel) for three bins in velocity dispersion as 
indicated by the labels. Note that the units of environmental density provide a relative measure of local density and are not associated to astrophysical objects 
such as field, groups or clusters. Right-hand panel: Fraction of rejuvenated galaxies as a function of galaxy velocity dispersion and dynamical mass (see top 
x-axis) for three different environmental density intervals —5 log p < — 1 (blue line), — 1 ^ log p < (orange line), and log p < 1 (magenta line). 
Error bars are Poisson errors. The figure shows clearly that the rejuvenation fraction depends on both galaxy mass and environmental density. 



environmental densities, while the intrinsic spread in both metal- 
licity and a/Fe is virtually independent of the environment. 

The slopes and zero-points for the scaling relations found in 
this work are in reasonable agreement with other recent results in 
the literature (see references in the Introduction). The slope of the 
age-cr relation is somewhat steeper than what we have found in 
T05. In particular the results for the high-density sample in T05 
are different, which might well be caused by a particularity of the 
Coma Cluster that dominates the high-density part of the sample 
in T05. Still, the impact on the derived formation epochs and their 
dependence on galaxy mass is minor (see Section|4j- The most im- 
portant difference with T05 and other previous work is that the stel- 
lar population scaling relations are independent of the environment 
for the bulk of the population. 



4 EPOCHS OF EARLY-TYPE GALAXY FORMATION 

The luminosity-weighted ages and a/Fe element ratios derived 
here enable us to approximate star formation histories. This is 
possible as the element ratio constrai ns formation time-sc ales. By 
means of a chemical evolution model dThomas et ai] 19991) . in T05 
we determined a simple scaling between a/Fe ratio and formation 
time-scales. The latter together with the average age defines a star 
formation history. The resulting star formation rates as functions of 
lookback time for various mass bins are shown in Fig. [9] Masses 
are derived from velocit y dispersion and effecti ve radius through 
the scaling provided by Cappel lari et al.l {2006). As discussed in 
detail in T05, these star formation histories are meant to sketch the 
typical formation history averaged over the entire early-type galaxy 
population (at a given mass). Real star formation histories of indi- 
vidual objects are expected to be more bursty and irregular. 

4.1 Comparison with T05 

The formation epochs of early-type galaxies below ~ 10 11 Mq 
are in good agreement with the values derived in T05 for the low- 
density sample. At the highest masses the situation is more com- 



plex. The steeper slope of the age-a relationship found here with 
respect to the low-density sample of T05 leads to significantly 
larger formation ages of the most massive early-type galaxies with 
log a ~ 2.5 (~ 10 12 Mq) by about 2 Gyr. The lookback time for 
their formation is ~ 12 Gyr (instead of ~ 9.5 Gyr in T05) and 
hence conspires to be very similar to the ages of the high-density 
massive early- types in T05. Their formation redshifts, however, are 
lower (~ 4) than in T05 (~ 5) because of the lower values for both 
matter density and Hubble parameter adopte d in the present work 
following the mos t recent determinations dTegmark et alj 120061 
IPercival et al]|2007h . 

4.2 Archaeological downsizing 

Fig. [9] illustrates the key observation that more massive galaxies 
form at higher redshift, now commonly referred to as 'do wnsiz- 
ing'. This naming was first suggested bv lCowie et alj d 19961) who 
studied deep survey data and found that 'the maximum rest-frame 
K luminosity of galaxies undergoing rapid star formation has been 
declining smoothly with decreasing redshift' . The present approach 
based on the 'archaeology of stellar populations', instead, mani- 
fests downsizing for the local galaxy p opulation, which logic ally is 
called 'archaeological downsizing' by Nei stein et alj JioOq) . Note 
that the latter is a universal ch aracteristic of galaxy formation and 
is not restricted to early-types {Gavazzi et alj|2002l; iHeavens et al.l 
12004 ICid Fernandes et alJr2007l : IPanter et alj|2007l) . Most impor- 
tantly, the convergence of local and high redshift observations to 
the general downsizing pattern is vital to set stringent constraints 
on the theory of galaxy formation. 

Clearly, a mechanism is needed to regulate star formation 
in massive objects in order to produce the anti-hierarchical trend 
suggested by observations. A GN feedback has been spotted as 
the most promisin g solution dSilk & Reesl [19981; lBinnev| |2004|; 
Granato et al .120041; Tsilkl2005l; iMiller et alj2006l;ISchawinski et al.1 
20061; iHopkins et alj I2OO8I ; iRenzinil I2OO9I ; ICattaneo et alj 120091 
Kormendv et al.1 120091). In a companion paper exploiting the 
MOSES sample dSchawinski et alj2007bl ; ISchawinski et al.l2 009a) 
we have found first empirical evidence that such a mechanism may 
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Figure 9. Specific star formation rate as function of look-back time for early-type galaxies of various masses as indicated by the labels. The grey hatched curves 
indicate the range of possible variation in the formation time-scales that are allowed within the intrinsic scatter of the a/Fe ratios derived. No dependence on 
environmental density is found. The upper x-axis connects time and redshift adopting Q m = 0.24, = 0.76, and Ho = 73 km/s/Mpc. Note that these star 
formation histories are meant to sketch the typical formation history averaged over the entire galaxy population (at a given mass). Real star formation histories 
of individual objects are expected to be more bursty and irregular. Intermediate- and low-mass galaxies in low-density environments get rejuvenated via minor 
star formation events below redshift z ~ 0.2 (see Section 3.8). This suggests a phase transition from a self-regulated formation phase without environmental 
dependence to a rejuvenation phase, in which the environment plays a decisive role possibly through galaxy mergers and interactions. 



indeed be occurring in early-type galaxies also at recent epochs. 
Moreover, AGN outflows at high redshift provide observational ev- 
idence for AGN feedback at early epochs in the evolution of galax- 
ies (e.g. lNesvadba et al.l2 008). The most recent renditions of semi- 
analytic, hierarchical ga laxy formation model s include a prescrip- 
tion of AGN feedback jPe Lucia et al . 2006; Croto n et al"]|2006l ; 
Bower et alj2006l;|Monaco et al.ll2007l;ICattaneo et alj|2005l.l2006l: 
Dekel & Bimboimll2006l; ICattaneo et alj|2007l: ICavaliere & Mend 
20071) . |Pe Lucia et al.l feOOot) show that the simulations now pro- 
duce star formation histories that are much closer to the obser- 
vational constraint presented in T05 (and Fig. [9] of the present 
work) than previous generations of semi-analytic models. The ob- 
served a/Fe enhancement of early-type galaxies and in particu- 
lar its correlation with galaxy mass (Fig. [3} needs yet to be recon- 
ciled with hierarchical models of galaxy fo rm ation dThomasll 19991 ; 
iNagashima etai]|2005l ; IPipino et al.ll2009l) . In IPipino et all d2009l) 
it is shown that the predicted a/Fe-mass relationship is still flat- 
ter and has considerably more scatter than the observations even in 
models with A GN feedback. This co nclusion is confirmed by the 
recent study of lArrigoni et alj ( I2009I) . who show that the a/Fe-o 
relation can only be produced if a top-heavy initial mass function 
and a lower fraction of binaries that explode as Type la supernovae 
are adopted. 



4.3 Phase transitions in the local universe 

A fundamental difference with respect to previous results is the lack 
of a dependence on environment. Fig. 10 in T05 presenting the for- 
mation epochs shows a significant delay in the formation of mas- 
sive galaxies in the field. As discussed in T05, this suggested that 
massive galaxies in low densities form later but on the same time- 
scales as their counterparts in clusters (see Fig. 10 in T05) as the 
direct consequence of early-type galaxies in denser environments 
having higher ages but the same a/Fe ratios. This implied that star 
formation must have essentially been on hold for the first 3 Gyr 
in massive objects in low densities. The new (and statistically more 
robust) results presented here suggest a different picture. In the new 
version of the star formation history plot (see Fig. [9) the main for- 
mation epochs of early-type galaxies as a function of galaxy mass 
are independent of the environmental density. 

However, some early-type galaxies are rejuvenated, i.e. they 
must have harboured minor star formation events at recent epochs 
within the past few Gyrs (see Section 3). Such recent star forma- 
tion activity that we call 'rejuvenation' occurs on top of the star for- 
mation histories shown in Fig. [9] The fraction of this rejuvenated 
galaxy population increases both with decreasing galaxy mass and 
with decreasing environmental density (see Fig. [8). This implies 
that t he impact of envi ronment increases with decreasing galaxy 
mass dTasca et alj2009l) . The dependence on environment suggests 
that galaxy interactions and merger activity might have been the 
major triggers of early-type galaxy rejuvenation in the past few 
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billion years since approximately redshift z ~ 0.2. During most 
of the Hubble time, instead, the formation and evolution of early- 
type galaxies was driven mainly by self-regulation processes and 
intrinsic galaxy properties such as mass. The stellar populations 
in early-type galaxies did not seem to feel in which environment 
the galaxy was forming. It is during this formation phase before 
z ~ 0.2 when the bulk of the stellar populations in early-type 
galaxies formed, and when the stellar population scaling relations 
(see Eqns.Q]and[2]( were established. This stage is then followed by 
a rejuvenation phase at redshifts below z ~ 0.2 affecting mostly 
low- and intermediate mass early-type galaxies, during which the 
environment must have played a key role possibly via galaxy merg- 
ers and interactions. 

Thus, there must have been a phase transition a few bil- 
lion years ago in the evolution of early-type galaxies from an 
environment-independent self-regulated phase of formation to an 
environment-dependent phase of (possibly) merger-driven rejuve- 
nation as sketched in Fig. [9] The universe may eventually move 
into a secular evolution phase in which intrinsic galaxy properties 
re-gain dominanc e over galaxy interactions as the main driver of 
galaxy evolution ( Kormendv & Kennicutt 2004). 



5 DISCUSSION 

In this paper we investigate the effect of environmental density on 
the formation epochs of early-type galaxies. The major improve- 
ment with respect to T05 and other previous work in the literature 
is the significant increase in sample size allowing for a proper sta- 
tistical analysis. We study 3,360 elliptical galaxies in a very narrow 
redshift bin (0.05 ^ z ^ 0.06) that have been morphologically 
selected from the SDSS DR4. The most important result is that, at 
fixed galaxy mass, we find no dependence of the stellar population 
parameters on environmental density for the bulk of the population. 
It is the fraction of rejuvenated ellipticals (about 10 per cent in the 
average environment), instead, that increases strongly with decreas- 
ing environmental density. In the following we compare this finding 
to other recent work in the literature. 

As outlined in the Introduction, pre-SDSS studies of the stel- 
lar population parameters in early-type galaxies based on rela- 
tively small, local samples consistently find younger average ages 
in low density environments. T05 have found an average shift in the 
luminosity-weighted age by ~ 2 Gyr. The present sample based 
on SDSS is significantly larger and therefore less biased towards 
the local volume, which may explain the discrepancy with previous 
work. The increase in sample size further implies a higher statisti- 
cal significance of the results. 

A first step toward lar ge samples b e yond the very local uni- 
verse has been done by iBernardi et al.l j 19981) who analyse the 
Mg 2 -<r relationship of a sample of 931 early-type galaxies in var- 
ious environments. Interestingly, they find the influence of the 
environment to be negligible, in agreement with the present re- 
sults. A more detaile d comparison is not possible, however, as 
IBernardi et afl GUI) 

do not derive stellar population parame- 
ters directly. In a mo re recent study based on the SDSS, instead, 
IBernardi et al. (2006) discuss the stellar population scaling rela- 
tions as a function of environment. Their Fig. 14 clearly agrees 
with our conclusion that these relationships are insensitive to the 
enviro nmental density. 

ISmith et alj (2006) analyse 3,000 red-sequence galaxies in 
nearby clusters and find some variation of the stellar population pa- 
rameters along the cluster radius. Early-type galaxies appear about 



0.06 dex younger and slightly less a/Fe enhanced in the outskirts. 
This result is not necessarily inconsistent with the present study. 
This effect on the average quantities is small and might well be 
driven by the increase in the relative number of rejuvenated ellipti- 
cals rather than by a change in properties of the bulk population. 

Interesting in the context of this paper is also the work by 
lHaines et al] d2006h based on SDSS DR4. In essence they find that 
the typical ages of dwarf galaxies show strong dependence on en- 
vironment, from being relatively old, passive satellites in clusters 
to being actively star forming and rejuvenated systems in the field. 
This is very well in line with our finding that it is the rejuvena- 
tion fraction (which is largest in the smallest galaxies) that strongly 
depends on environment. 

Colour gradients (with bluer colours at larger radii) are 
found to become steeper with decreasi ng environmental density 
dLaBarberaetal.ll2005l : iKo & Imll20o3) . Combined with the fact 
that colour gradients are produced by meta l licity rather than by 
age JSaglia et al.l l200d : iTamura et alj|2000l : iMehlert et al] 120031 : 
IWu et al. l2005h . this suggests that a higher (dry) merger rate in 
denser environments might have flattened existing metallicity gra- 
dients. However, the results of the present work imply that the for- 
mation and evolution of the stellar populations in galaxies must still 
be controlled intrinsically independent of the environment. Such 
dissipation-les s, dry mergers in dense environments ha ve indeed 
been observed ( Ivan Dokkum et al.ll999l ; lTran et alj2005l) . But they 
must occur relatively late in the evolutionary history of the early- 
type galaxy population and must be moderate in number in order 
to keep the stellar population scaling relations intact. The obser- 
vationally constrained, modest merger rate for massive galaxies of 
one merger per galaxy since z ~ 1 may be con sistent with that 
requirement (Bell et al. 2006; Pozzetti e tai]|20ld) . 

Most recently, Cooper et al.l ([2009) have found a residual age- 
density relationship for red sequence galaxies drawn from SDSS, 
which is in apparent contradiction to the results of this paper and 
other SDSS based stellar population studies such as lBernardi et al.l 
(2006). These two results can be reconciled, however. It is possible 
that the actual variation in age as a function of environment is too 
small for a direct detection in the prese nt work, but visible as the 
residual age-density relation discussed in ICooperetaljj2009l) . This 
would be consistent with the work of Clemens et al (2006, 2009), 
who find a small (< 0.1 dex) systematic offset toward older ages 
in high density ellipticals independently of galaxy mass. 



6 CONCLUSIONS 

We analyse the stellar population properties of 3,360 early-type 
galaxies as functions of their environmental density and intrinsic 
parameters such as velocity dispersion and dynamical mass. A ma- 
jor aim is to check the result found in T05 by increasing the sample 
size by more than a factor ten. The objects are morphologically se- 
lected by visual inspection from the SDSS in the narrow redshift 
range 0.05 ^ z ^ 0.06. The sample utilised here is part of a 
project called MOSES: Morphologically Selected Early-types in 
SDSS (see also lSchawinski et alj|2007br) . The most radical differ- 
ence with respect to other galaxy samples constructed from SDSS 
is our choice of purely morphological selection. We have visually 
inspected 48,023 objects in the redshift range 0.05 ^ z ^ 0.1 with 
r < 16.8 mag and divided the sample in 31,521 late-type (spi- 
ral arms, clear disc-like structures) and 16,502 early-type galax- 
ies. The major advantage of this strategy is that our sample is not 
biased against star forming early-type galaxies. We calculate a 3- 
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dimensional volume density at the location of each object, provid- 
ing us with a local environmental den sity. Kinematics o f gas and 
stars are determined using th e code by Sarzi et al. (2006) which is 
based on the ppxf method bv lCappellari & Emsel lem (2004). Stel- 
lar population and emission line templates are fitted simultaneously 
to the galaxy spectrum, which yields clean absorption line spec- 
tra that are free from emission line contamination. Subsequently, 
Lick absorption line i ndices are measured. By m eans of the stellar 
population models o f lThomas et ail J200il2004l) the stellar popu- 
lation parameters age, metallicity, and a/Fe ratio are derived using 
a minimised \ 2 technique. Dynamical masses are calculated from 
a scaling between gala xy velocity dispersion and half-light radius 
dCappellari etai]|2006h . 

The morphological selection algorithm is critical as it does 
not bias against recent star formation. We find that the distribution 
of ages is bimodal with a major peak at old ages and a secondary 
peak at young ages around ~ 2.5 Gyr containing about 10 per cent 
of the objects. This is analogue to 'red sequence' and 'blue cloud' 
identified in galaxy populations usually containing both early and 
late type galaxies. The fraction of the young, rejuvenated early-type 
galaxy population depends on intrinsic galaxy properties as well as 
their environmental densities, i.e. it increases as a function of both 
decreasing galaxy mass and decreasing environmental density. It 
varies from only a few per cent for the most massive galaxies in 
densest regions up to 45 per cent for the smallest early-type galax- 
ies (a few 10 10 Mq) in the loosest environments. This implies that 
the impact of environment increases with decreasing galaxy mass. 

The rejuvenation process detected here through optical ab- 
sorption line indices may date back or stretch over several Gyrs 
back in time. Most likely an early-type galaxy has been rejuvenated 
through several minor star formation events each involving star for- 
mation of only a few per cent of the galaxy's mass. Most impor- 
tantly, the rejuvenated galaxies display slightly lower a/Fe ratios 
and slightly higher total metallicities. The former supports the reju- 
venation scenario, as a depleted a/Fe ratio is a clear indication for 
late Type la supernova enrichment caused by residual star forma- 
tion events. This is further supported by the fact that the emission 
line spectra of most of the rejuvenated galaxies indicate the pres- 
ence of minor, ongoing star formation. Some fraction of these are 
star forming/AGN composites. All objects that host AGN in their 
centres without any star formation being involved have stellar old 
populations. Finally, the higher metallicities suggest that the latter 
does not involve purely metal-poor pristine gas but must include 
at least some fraction of previously enriched interstellar medium, 
either from internal re-processing or accretion. 

The early-type galaxy population on the red sequence, in- 
stead, turns out to be insensitive to environmental densities. We 
find that the stellar population scaling relations including the age- 
mass relation are independent of the environment and only driven 
by galaxy mass. We confirm and statistically strengthen earlier re- 
sults that luminosity weighted ages, metallicities, and a/Fe ele- 
ment ratios correlate with velocity dispersion. We further show that 
these 'stellar population scaling relations' are well defined also for 
galaxy mass, which reinforces the now widely accepted conception 
of 'archaeological downsizing' of early-type galaxies. Most inter- 
estingly, however, these scaling relations are not sensitive to envi- 
ronmental densities and only driven by galaxy mass. This and the 
old average ages of early-type galaxies suggest that the formation 
of early-type galaxies is environment independent and driven only 
by self-regulation processes and intrinsic galaxy properties such as 
mass. This formation phase has lasted most of the Hubble time, 
but appears to have now come to an end. At lookback times of a 



few billion years, below z ~ 0.2, galaxy interactions and (minor) 
mergers appear to have become the major players in galaxy evolu- 
tion leading to rejuvenation of mostly intermediate and low-mass 
early-type galaxies in low density environments. 

We infer that early-type galaxy formation has undergone a 
phase transition a few billion years ago around z ~ 0.2. A self- 
regulated formation phase without environmental dependence has 
recently been superseded by a rejuvenation phase, in which the en- 
vironment plays a decisive role possibly through galaxy mergers 
and interactions. The universe will subsequently move into a sec- 
ular evolution phase in which intrinsic galaxy properties re-gain 
dominanc e over galaxy interactions as the main driver of galaxy 
evolution dKormendv & Kennicuttll2004l) . 

The galaxy catalogue produced in this study can be found at 
w w w. icg . port . ac . uk/~ thomasd. 
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